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Abstract. This work addresses the problem of estimating the parameters of the general 
half-normal distribution. Namely, the problem of determining the minimum risk equiva¬ 
riant (MRE) estimators of the parameters is explored. Simulation studies are realized 
to compare the behavior of these estimators with maximum likelihood and unbiased es¬ 
timators. A natural Monte Carlo method to compute conditional expectations is used 
to approximate the MRE estimation of the location parameter because its expression in¬ 
volves two conditional expectations not easily computables. The used Monte Carlo method 
is justified by a theorem of Besicovitch on differentiation of measures, and has been slightly 
modified to solve a sort of “curse of dimensionality” problem appearing in the estimation 
of this parameter. This method has been implicitly used in the last years in the context 
of ABC (approximate Bayesian computation) methods. 
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1 Introduction 


Let Z be a A/'(0,1) random variable. The distribution of X := \Z\ is the so-called half-normal 
distribution. It will be denoted iLA/'(0,1) and its density function is 



A general half-normal distribution HN{^^r]) is obtained from HN{0^ 1 ) by a location-scale trans¬ 
formation: HN{^, rf) is the distribution of F = ^ + tjX. 

The classical paper Daniel (1959) introduces half-normal plots and the half-normal distribution, 
a special case of the folded and truncated normal distributions (see Johnson, Kotz and Balakrishnan 
(1994)). Bland and Altman (1999) and Bland (2005) propose a so-called half-normal method to 
deal with relationships between measurement error and magnitude, with applications in medicine. 
Pewsey ( 2002 ) uses the maximum likelihood principle to estimate the parameters, and presents a brief 
survey on the general half-normal distribution, its relations with other well-known distributions and 
its usefulness in the analysis of highly skew data. Pewsey (2004) proposes bias-corrected versions of 
the maximum likelihood estimators. Nogales and Perez (2011) deals with the problem of unbiased 
estimation for the general half-normal distribution. 

Here we consider the problem of equivariant estimation of the location and scale parameters, ^ 
and 7 ^, but first we provide a brief review of results for unbiased and maximum likelihood estimation 
appearing in the literature. 

The density function of HN{^,r]) is 


fviy) = -fx 
V 


1 ( y-^ 

2 V V 


It is readily shown that 


/ 2 jr _ 2 

E(Y) = ^^ri\— and Var(F) = - 77 ^. 

V TT TT 

Let us recall a lemma from (Nogales and Perez (2011)). We write for the standard normal 
cumulative distribution function. 

Lemma 1. Let = |Z^|, 1 < i < n, where Zi,..., is a sample of the standard normal distribution 

X(0,1). Let Cn := X(Xi:^), where Xi;^ denotes the minimun of Xi,..., X^. 

(i) Cn = /o°°(2 - 2$(t))” dt. 

(ii) For n > 1, c„ < + ^). 

Let Fi,..., W be a sample of size n from a general half-normal distribution with unknown parame¬ 
ters, ^ and 77 . Yi-n denotes the minimum of Yi,..., F^. From the factorization criterion, we obtain that 
{T.UY^,T.UY„Y,..n) is a sufficient statistic. Indeed, it is minimal sufficient, although not com¬ 
plete. With the notations of the lemma, we write Yi = ^-h 77 X^. Notice that Fi;^ = min^ Yi = ^ + 77 X 1 :^ 
and E{Yi,n) = i Y r]Cn. 

The next proposition (Nogales and Perez (2011)) yields unbiased estimators of the location and 
scale parameters, ^ and 77 . Both estimators are L-statistics and functions of the cited minimal sufficient 
statistic. 


Proposition 1. Let Fi,... ,F^ be a sample of size n from a general half-normal distribution with 
unknown parameters, ^ and 77 . 

(i) ^ := ig unbiased estimator of the location parameter 

V ^ 

(ii) 77 := unbiased estimator of the scale parameter 77 whose distribution does not 


depend on 
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Remark, We also have that the sample mean Y is an unbiased estimator of the mean ^ + 
Moreover, an unbiased estimator of rf is 



where 5^ ;= ^ Er=i(^* - ^ f is the sample variance; notice that its distribution does not depend 
on Y and S‘^ also are functions of the sufficient statistic given above. The reader is referred to 
Nogales and Perez (2011) for these and other results about unbiased estimation of the parameters of 
the general half-normal distribution. □ 


Remark, Pewsey (2002) provides maximum likelihood estimates for each of the parameters ^ and r]: 

/ ^ n \ 




i=l 


A large sample based bias-correction is used in Pewsey (2004) to improve the performance of the 
maximum likelihood estimators ^ and fj. □ 


2 A Monte Carlo method to approximate conditional expec¬ 
tations 


In this section, we describe a natural Monte Carlo method to compute conditional expectations based 
on a theorem of Besicovitch on differentiation of measures. It will be used in the next section to 
approximate the minimum risk equivariant (MRE) estimator of the location parameter ^ because its 
expression involves two conditional expectations not easy to compute. 

We first recall briefly a theorem of Besicovitch (1945, 1946) for differentiation of measures (see, 
for instance. Corollary 2.14 of Mattila (1995)). This theorem extend to Radon measures the classical 
Lebesgue Differentiation Theorem. 

Theorem 1 (Besicovitch (1945, 1946)). Let X be a Radon measure on and f :W^ a loeally 
X-integrable funetion. Then 

lim . . ^. .. [ fdX = f(x) 

NO X{Br{x)) Jb^{x) 

for X-almost all x where Br{x) denotes the ball of eenter x and radius r > 0 for the norm || • ||oo 

on 


Now let (11, M, P) be a probability space, X : (H, M, P) RP be an n-dimensional random variable 
and Y : (11,M, P) ^ M be a real random variable with finite mean. The conditional expectation 
E{Y\X) is defined as a random variable on R^ such that fx-i(B)^ ~ 

Borel set B in where P^ denotes the probability distribution of X. 

Although the existence of the conditional expectation is guaranteed via the Radon-Nikodym theo¬ 
rem, its computation is, generally, involved. Nevertheless, according to the previous result, for P^- 
almost every x G 




L 


Y (cd) dPiuj) = lim 

(Be(a=)) ^-^0 


L 


B,{,x) 


E{Y\X 


= x')dP^{x')=E{Y\X =x) 


By the Strong Law of Large Numbers, for almost every sequence (w,) in Q,, we have 




^ ^ l- 

and 


k 

[ E{Y\X = x')dP^{x') = 
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where I a denotes the indicator function of A. Observe that, for every e > 0, the rate of convergence 
is ll\/n. 

Hence, we have proved the following result: 

Theorem 2. Let be a probability spaee, X : {Vt^A^P) BP be an n-dimensional random 

variable and Y : (ft^A^P) B be a real random variable with finite mean. Then, for P^-almost 
every x G B^ and almost every sequenee (cji) in Ll, we have 


E{Y\X = x) 


= limlim 

e-j^O k 




This theorem yields a means of approximating the conditional expectation of Y given X. The 
following simple example illustrates the method. 

Example 1. Let (X, T) be a bivariate normal random variable with null mean such that Var (X) = 
Var (T) = 1 and Cov(X, T) = 0.5. In this case, there is no need for an approximation to the 
conditional expectation of Y given X = x because it is x/2. The conditional distribution of Y given 
X = X is X(|x, |\/3)- Applying the proposed method to evaluate E{Y\X = 1), given a small 
e > 0, we may choose a sample {xi,yi)i<i<k from the joint distribution of X and Y and approximate 
E{Y\X = 1) by 


Taking e = 0.1,0.01 and samples from the joint distribution of X and Y with sample sizes k large 
enough to obtain m = m{k) = Yli=i ^[i-e,i+e] (^i) = IOO 5 1000, 5000, we obtained the approximations 
for E(Y\X = 1) summarized in Table 1 and Figure 1; 100 replications of each simulation have been 
conducted to obtain the table and the figure. Namely, taking m = 1000, for instance, the value 
0.493947 appearing in the table as an approximation of E{Y\X = 1) when e = 0.1 is the mean of the 
100 values of the quotient (1) obtained after 100 replications of the experiment of choosing a k-sized 
sample {xi,yi)i<i<k of the joint distribution of (X, T), k being large enough to get m = m{k) = 1000. 
Table 1 also includes the “mean squared error” (MSE) calculated from these 100 values: the format 
used for a typical entry in the table is E{Y\X = 1) ± MSE. The box-plot of the figure describes the 
distribution of these 100 values (a dotted red line represents the mean). 


m 

100 

1000 

5000 

e = 0.1 

0.505885 & 0.006395 

0.493947 & 0.000815 

0.497892 & 0.000128 

e = 0.01 

0.503655 & 0.007165 

0.499826 & 0.000716 

0.499471 & 0.000150 


Table 1. Approximation of E{Y\X = 1) & MSE as a function of the number of simulations, m, for 
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Remark, The described method of Monte Carlo approximation to the conditional expectation E{Y\X = 
x) is based on the naive idea that one can approximate it from a sample {xi,yi)i<i<ri by the mean 
of the Hi corresponding to points Xi lying in a narrow neighborhood of x. From a probabilistic point 
of view, the method has been justified by the mentioned theorem of Besicovitch on differentiation of 
measures. When the joint density of X and Y is known, E(Y\X = x) is the mean of the conditional 
distribution of Y given X = x, and the problem of compute a conditional expectation is reduced to 
the problem of computing a mean. Notice that the existence of a joint density is not required by the 
method and it could be specially useful when densities are not available or are not easy to compute 
(see the next example). □ 

Example 2. (Example 1, continuation) A similar simulation study has been performed to approxi¬ 
mate the conditional expectation E{y\U = 0.5), where V = sin(X • Y) and U = cos(X^ + F^); the 
obtained results are: 


m 

100 

1000 

5000 

6 = 0.1 

0.127650 & 0.001890 

0.127280 & 0.000202 

0.124169 & 0.000025 

6 = 0.01 

0.123063 & 0.001620 

0.125869 & 0.000153 

0.1252856 & 0.000031 


Table 2. Approximation of E{V\U = 0.5) & is the sample variance) as a 

function of the number of simulations, m, for e = 0.1,0.01. 



Figure 2. Box plots of the approximations of E{V\U = 0.5) 
as a function of the number of simulations, m, for e = 0.1 and e = 0.01. 

Remark, In a classical statistical framework, we can provide additional guarantees on the method, 
since the obtained Monte Carlo approximation to the conditional expectation E{Y\X = x) coincides 
with the value at the point x of the kernel estimator (the Nadar aya-Wat son estimator) of the regression 
curve y = E{Y\X = x) for the kernel K{x) = /[_i^i](x) (see Nadaraya (1989), p. 115). From this 
point of view, e plays the role of the bandwidth parameter. We refer to Hardle (1992, Ch. 5) for a 
detailed discussion on the important problem of the choice of the bandwidth. □ 

Remark, As it is pointed out to us by the referees, in a Bayesian setting a similar idea has been in use 
in recent years to generate an approximate sample from the posterior distribution given x assuming 
that the likelihood function is easy to sample. This proceeds by sampling values Oi from the prior 
distribution and Xi from the distribution of the data given and accepting those parameters Oi such 
that Xi is in the ball B^{x) centered at x of radius e > 0. In fact, from a sample of size k we can 
approximate the posterior probability given x of a subset T of the parameter space by 

Ya= 1 lB,(x){Xi)lT{0i) 

We also can approximate the posterior mean given x of a function / of the parameter by 
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Remark, In this paper, the main application of the Monte Carlo method for the approximation of 
conditional expectations is given in the next section to approximate the estimation of the location 
parameter of the general half-normal distribution, because it is defined in terms of a quotient of two 
not-easily-computable parameter-free conditional expectations given a (n — l)-dimensional statistic U. 
Some “curse of dimensionality problem” appears when n is large because, in this case, it is not easy 
to find large samples of points lying in a small ball centered at a point U{y). This is why we had to 
modify the Monte Carlo method for the approximation of conditional expectations taking advantage 
of the underlying distribution of Y (the general half-normal distribution) and the invariance properties 
of U. This could become an important scholium of the paper, as the ideas used here could be useful 
to deal with the “curse of dimensionality problem” in similar situations. □ 

3 Equivariant estimation of the location parameter of the ge¬ 
neral half-normal distribution 

In this section we consider the problem of determining the minimum risk equivariant estimator of the 
location parameter ^ of the general half-normal distribution HN{^^r]) when the scale parameter 77 is 
unknown. We cannot provide an explicit expression for this estimator, since it is described in terms 
of two conditional expectations that had to be approximated by simulation. 

To achieve this goal, an R program was developed based on the method of computing conditional 
expectations described in the previous section. In fact, the method has been slightly modified to solve 
a sort of “curse of dimensionality” problem. 

We consider the scale-location family of densities 

f \ ^ yn-C\ 


where 

f{yi,-,yn) = 0^ expj-i^y^l 

This family remains invariant under transformations of the form ga,b{yi^ •••^yn) = {ciYhyi^ ..., aYhy^)^ 
a G M, 6 > 0. 

To estimate the location parameter ^ when the scale parameter 77 is unknown, we have the next 
result, a direct consequence of classical equivariant estimation theory (see Lehmann (1983)). First, 
recall that an estimator T of the location parameter is equivariant if T(a + 6 x 1 ,..., a + hxn) = 
a + 6 T(xi,... ,x^), for all a G M and all b > 0. 

Proposition 2. When the loss funetion W 2 {x; = r]~^{x — ^)^ is eonsidered, the MRE estimator 

i ofi is 

e = T*-(pot/)T* 

where 


n 

T^ = Y, = 

n . 

1=1 

( Y^-Yn Yn-2 - Yn F„-l - \ 

'vFn-i - ’ ■ • ■ ’ ’ |r„_i - r „|7 ’ 

i;5=o,,=i(Ty|c/) 

Remark, Tq can be replaced by any other equivariant estimator of ^ (i.e., satisfying To*(a+ 67 /i,..., a+ 
byi) = a -h 6 Tq (7/1,... ,7/1) for every a G M, 6 > 0 ), and can be replaced by any positive estimator 
of 77 satisfying Tf (a + 67/1,..., a + byi) = bT^ (7/1,..., 7/1) for every a G M, 6 > 0 . □ 
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A simulation study has been performed to investigate the behavior of the minimum risk equivariant 
estimator In it, we used 100 simulations with sample sizes n = 100,1000, 5000 from the HN{W^A) 
distribution, obtaining the results summarized in Table 3 and Figure 3 (see below how we have made 
use of the method of approximation of conditional expectations to obtain the values of the Tables 3 
and 4). 


m 

100 

1000 

5000 

e = 0.1 

9.412753 & 2.307873 

9.517691 & 1.881063 

9.732626 & 0.158635 

e = 0.01 

9.603969 & 0.687243 

9.274164 & 3.501158 

9.867600 & 0.027826 


Table 3. Approximations of ^ & MSE as a function of the number of simulations, m, for e = 0.1,0.01. 




Figure 3. Box plots of the approximations of ^ & MSE 
as a function of the number of simulations, m, for e = 0.1 and e = 0.01. 

To compare the behavior of the unbiased estimator the maximum likelihood estimator ^ and the 
minimum risk equivariant estimator we used 100 simulations with sample sizes n = 100,1000, 5000 
from the HN{W^A) distribution, obtaining the results summarized in Table 4 and Figure 4: 



m 

100 

1000 

5000 

1 

e = 0.1 

9.997350 & 0.003289 

9.999356 & 0.000020 

10.000823 & 0.000001 

e = 0.01 

9.999662 & 0.003443 

9.999989 & 0.000025 

10.001050 & 0.000002 


e = 0.1 

10.047256 & 0.005455 

10.004383 & 0.000039 

10.000823 & 0.000001 

e = 0.01 

10.049128 & 0.005752 

10.005005 & 0.000050 

10.001050 & 0.000002 


e = 0.1 

9.412753 & 2.307873 

9.517691 & 1.881063 

9.732626 & 0.158635 

e = 0.01 

9.603969 & 0.687243 

9.274164 & 3.501158 

9.867600 & 0.027826 


Table 4. Approximations of ^ & MSE, ^ & MSE and ^ & MSE as a function 
of the number of simulations, m, for e = 0.1, 0.01. 


Table 4 and Figure 4 illustrate the biased character of the maximum likelihood estimator ^ and 
the minimum risk equivariant estimator 

Let us describe in more details the ideas used in these simulations. For a sample y = (^i,..., 
n = 100,1000, 5000, of the distribution iLA^(10,4), we have 

N 

Kf/(y)) = lim^ 
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m = 100 






Figure 4. Box plots of the approximations of ^ and ^ as a function of the 
number of simulations, m = 100,1000, 5000, for e = 0.1 and e = 0.01. 
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where 


N,= f f{y')dy', D,= [ giy')dy\ 

J A^{y) J A^{y) 

/(y')=To*(y')Ta*(y')exp|-i||y'||^|, 9{y') =n{y'f e^v\-\\W\\l^ , 


Ae{y) = {y' e [O, lO]": max \Ui{y') - Ui{y)\ < e}. 
Now, take a sample S of A^{y) and approximate and by 

^ ^ y'es ’ y'es 


respectively. So, p{U{y)) can be approximated by 


Ciy) := 


T^v’esfjy') 

T^y'esyiv') 


and C(y) is approximated by D{y) := T^{y) - C{y)T^{y). 

To approximate C{y), a first idea would be to divide the interval [0,10] in multiple subintervals of 
small length e > 0 and consider the grid in the interval [0,10]formed by the n-power set of the ends 
of these subintervals (we have restricted ourselves to the interval [0,10] because the functions f{y) and 
g{y) are almost null when one of the coordinates of the vector y is greater than 10). The sample S 
would then be formed by the grid nodes that are in A^. The main problem with this approach is that 
the size m of the sample S is very small: it becomes smaller as n increases, because of the so-called 
“curse of dimensionality” problem. In order to avoid this problem and obtain a sample size m large 
enough for S (given n, we take m = lOOn), we have used the following algorithm, a modification of 
the described Monte Carlo method to approximate conditional expectations that hinges on the use of 
the invariance of U under scale and location transformations. Namely: 


Step A. Let n G N and be ^ = (^i,... ,^n) a n-sized sample of the distribution iLA/'(10,4). For 1 <i < 
n — 2, let ai := and take 0 < e < min{0.1, mini <^<^_2 |ni||. 

Step A.l. At this stage we choose coordinatewise at random 100 • n vectors = {v^\ ... 

1 < j < lOOn, in such that maxi<^<^_i \Ui{v(J^)-Uiiy)\ < € as follows: 

A.1.1. Make j = 1. 

A.1.2. Take at random in [0,10] such that has the same sign as yn-i—yn- 

(So, the last coordinates of U{v^A^ and U{y) are the same). 

A. 1.3. For 1 < i < n — 2 take at random on the interval determined by Vn^ + ('^n-i ~ 
Vn^){ai - e) and ~ + e). (So \Ui{v^^^) - Ui{y)\ < e). 

A.1.4. Make j = j + 1 a go back to Step A.l until lOOn vectors = [v[^\ ..., Vn'^), 1 < j < 
lOOn are obtained. 

Step A.2. Since the vectors v^A = 1 < J < lOOn, do not lie necessarily in [0,10]^ (so 

neither in A^{y))^ we can make some random location-scale transformations to put them into 
[0, lO]’^. These transformations do not modify the required fact that maxi<^<^_i \Ui{v^A^ — 
Ui{y)\<e. 

A.2.1. If < 0 for some we define = v-\-vP ^ l<i<n, l<j< lOOn, where v is 
choosen at random between — mini<^<^^i<j<ioonand 1 — mini<^<^^i<j<ioon 
Otherwise, uP = vP, I <i <n, 1 < j < lOOn. 

A.2.2. Each vector u^A ig divided by maxi<^<^i4p^ and multiplied by a random number 
choosen in [0,10] to obtain the vector w^A ^ 
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A.2.3. Take S = : 1 < i < lOOn} and approximate C{y) by 


^10°;/(«;«)) 


and D{y) by T^iy) - C{y)T,*{y). 

Step B. Finally, following the process designed in Step A, we choose k := 100 random samples of 
size n from the FfA'(10,4) distribution and approximate the mean and the mean squared error 
of i by 

k k 

and ^ ^(D(yW) - 10)^, 

^ i=l ^ i = l 

respectively, and we construct a box-plot with the values 


Remark, Notice that both ^ and ^ are equivariant estimators of the location parameter So they 
have greater risk for the loss function ld ^2 than Hence, in the previous simulation study, the MSE 
of ^ should have been smaller than the MSE of ^ and That has not been the case because, for the 
MRE estimator, we have not real estimates of but approximations of these estimates obtained by 
a modification of the Monte Carlo method of computing the conditional expectations appearing as 
the numerator and denominator of a quotient. But this is a possible issue to approximate minimum 
risk estimations of a location parameter, and a possible way to avoid the “curse of dimensionality 
problem”. □ 

Remark, Although less interesting from the perspective of real applications, for completeness we now 
consider the problem of estimating the scale parameter ^ when the location parameter rj is known, say 
T] = r]o. In this case, the joint density of Yi,..., is 

1 r2 ^ f 1 ^ ^ 1 

•••,?/«) = - exp<^ 1 -fK.+oo[(l/l:n), 

Vo y ^ ^Vo 


where yi:n := min{^i,..., This family remains invariant under translations of the form ga{yi,... ,yn) 
{yi- a,...,yn- a). 

The equivariant estimator of minimum mean squared error of the location parameter ^ is 




V2^ ^ 




In fact, for the loss function ~ MRE estimator of the location parameter ^ 

is the Pitman estimator 


Ti{yi,... ,yn) 


ufo{yi - u,yn - u)du 
Mvi - U,yn - u)du 


Eor y G we write y for the mean of ^i,... After some algebraic manipulations, we obtain: 


ufo{yi - u, ...,yn - u)du = 




exp 


/. 


-^exp <1 -^{yi:n - yf \ + y\/^^ {^{yi-.n - y) 


{■ 




Vo 
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and 



and the statement follows easily from these expressions. □ 


4 Equivariant estimation of the scale parameter of the general 
half-normal distribution 


Unlike what happens with the location parameter for the scale parameter r] an explicit expression 
for the MRE estimator is obtained. 

Recall that an estimator T of the scale parameter r] is equivariant if T(a + ..., a + bxn) = 

bT{xi, ..., Xn), for all a G M and all b > 0. 


Proposition 3. When using the loss funetion Wi{x;^,r]) = r] ‘^{x — r])‘^, the MRE estimator fj of r] 
is 


f] = 




n+1^ ^n+1 


2 


^n+2 


/ n(n+] 

V n-1 


+ 1) Y-Yi:, 


, OO 


n(n+2) 
n-1 S 


S'. 


where tn denotes the Student ^s t-distribution with n degrees of freedom, is the sample varianee and 
r denotes Euler’s gamma funetion. 


Proof. The MRE estimator of the scale parameter p, when using the loss function Wi, is 

o. . ^ |+°° v^f{vy[, vy'^_.^)dv 

where /' is the joint density when = 1 of ¥( := Yi—Yn, 1 < i < n — 1, and y'- := yi—yn, 1 < i < n — 1. 
Notice that 


/ + 00 

f{yi+t,...,yn + t)dt 

-OO 


' —OO 

n 


t i=l 


I 


exp 




-yi:n 

OO 


V^iy-yi-.n) 


{-^it + yf]dt 
[-lu^]du. 


exp 


Hence, for /c G N, applying Eubini’s Theorem after a suitable change of variables in the inner 
integral, 


1 

\/n 



POO 

4(y) ;= / v’^f'{vy[,...,vy'„_^)dv 

Jo 

1 ;'“ exp |-i(n-1)^2 


exp 


_ 1 /2\2 r 

Vn Y) 


Vniy-yi-.n) 


Vniy-yi:n) 

Jk{t,y)dt. 


{-^“1 


dudv 
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where 


Jk{t,y) := j v*+^exp + (n - 1)S'^(2/))| dv 


2fe/2r 


if + (n- 1)52 (y))''^'' 


where, for t > ^/n{y — ^i:n), we have made the change of variables w = + (n — l)S‘^{y)). 

So, 


Finally 


□ 


(lY 2>^/2r ( 


\/n V 


TT 

n + fc 


2^r(^) 


(n — 1) 2 S{y)^+^ 


'k + 2\ dt 

) Js/^in.r.) (i2 + (n-l)52(y))^ 

^fc+i 


n{k + 1) y - yi:„ 
n-1 5(y) 


-, 00 


f]{y) = 


injy) 

In+liy) 


-ir(^) 


n±l\ ^n+1 


2 r(^i±2) 


tn+2 


nin-\-l) y-yi:n 

n-1 S(y) ’ 


OC 


n(n+2) y-yi:n ^ 
n-1 Siy) 


¥rS{y). 


Remark. A simulation study has been performed to compare the behavior of the unbiased estimator 
fj^ the maximum likelihood estimator fj and the MRE estimator f] using 1000 simulated random samples 
of size n = 10,20,30 from the iJA^(10,4) distribution. The results obtained for the means and the 
mean squared errors of the three estimators are presented in Table 5 and Figure 5 (as before, a dotted 
red line represents the mean). 


n 

fj 

fj 

V 

10 

3.996009 & 1.052443 

3.520680 & 0.952987 

3.568520 & 0.929288 

20 

3.996575 & 0.526328 

3.760888 & 0.458780 

3.795590 & 0.450882 

30 

4.015727 & 0.324161 

3.845478 & 0.294937 

3.871677 & 0.291209 


Table 5. Sample mean and MSE of the estimators calculated using 1000 random 
samples of size n from the HN{10,4) distribution. 



^10 20 ^ 
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Figure 5. Box plots for the estimator fj for sample sizes n — 10, 20, 30 (above) 
and for the estimators fj, f) y f) for sample sizes n = 10, 20, 30, respectively (below). 


Notice that both fj and fj are equivariant estimators of the scale parameter r]. So they have greater 
risk for the loss function Wi than fj. Hence (see Table 5 and Figure 5), in the previous simulation 
study, the MSE of fj is smaller than the MSE of fj and fj. □ 

Remark. Although less interesting from the perspective of real applications, for completeness we now 
consider the problem of estimating the scale parameter 77 when the location parameter ^ is known, say 
^ = ^0- After the shift ( 7 / 1 ,. .., 7/^) 1-^ (^1 — Co? • • • 5 ~ Co)? the statistical model remains invariant 

under the transformations (dilations) of the form ( 7 / 1 ,, 7/^) (<^^i? • • • ? <^^n)? for a > 0. Eor the loss 

function W{{r],x) = {x — rjY/ri‘^^ the MRE estimator of the scale parameter 77 is 


where B denotes Euler’s beta function. In fact, for the loss function the MRE estimator of 77 is 

poo 

/ v^hi{v{yi - ^o),...,v{yn - ^o))dv 
Jo 


^2(^1? •• • ?^n) — 


— 


poo ? 

/ - ^o),-,viyn -^o))dv 

Jo 


where 


hl{yi,...,yn) = (“) -2^yi p[0,+oo[(2/l:n)- 


i=l 

To simplify the notation, we assume without loss of generality that Co = 0- The change of variable 

b 


t = \ leads to, ioi k = n,n ^ 1, 


f 


v’^hi{vyi,... 


n-\-k— 1 

,vyn)dv = 2 2 TT” 



fc + 1 


^ [0, + OO 


[{yi-.n): 


and the assertion then follows easily. 

Note also that, when C = Co? 

i=l 

is the minimum variance unbiased estimator of 77^. This is a consequence of the Lehmann-Scheffe 
Theorem and the facts that “Co)^ is a sufficient and complete statistic and “Co)^ 
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has a x^{n) distribution. A little more work shows that 


r(f) 


i=l 


g(t4) 


- ^o)2 

i=l 


is the minimum variance unbiased estimator of 77. □ 
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